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<N 1 Abstract 

The known analytic properties of the Compton amplitude at small Q place significant constraints on 
Q | its behaviour at large Q 2 . This calls for a re-evaluation of the role of perturbative evolution in past 
fits to data. 
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There is a widespread belief that the rise at small x observed at HERA in the proton structure function 
is associated with the presence of a singularity at N = in the Mellin transform of the DGLAP splitting 
function. Indeed, fits to structure functions typically incorporate perturbative QCD evolution in which 
this singularity plays a key roleM ; but ignore an oldt 2 ! but very basic piece of knowledge: that at small 
x a structure function is governed by Regge theory! 3 ] . Recently, two of us have made^ an extremely 
successful fit to all small- x data, from Q 2 = up to 2000 GeV 2 , in which we concentrate on the Regge 
theory and ignore the requirements imposed by perturbative evolution. It is the purpose of this paper 
to try to bring the two together. In particular, we argue that, by analytically continuing in Q 2 , one 
can conclude that the singularities in the complex A-plane of the Mellin transform of the structure 
function must already be present at small Q 2 , and the perturbative evolution cannot generate new 
singularities that appear only at high Q 2 . In particular, the role customarily assigned to the N = 
singularity is not sustainable. 

Regge theory has its basis in the analyticity properties of the appropriate scattering amplitude, in 
this case the virtual Compton amplitude whose imaginary part at t = contains the familar structure 
functions Fi(x,Q 2 ) and F2(x,Q 2 ). A rigorous analysis of analyticity properties of scattering ampli- 
tudes is extremely difficult t 5 l So instead one assumes^ that the region of analyticity is the same as 
the analyticity region shared by all the infinite number of Feynman graphs that one can draw for the 
amplitude, ignoring the fact that in a strong-coupling regime the numerical values of these graphs are 
irrelevant because the perturbation series diverges. From this approach one may deduce that, both 
for the real-photon Compton amplitude and for Q 2 < 0, there is sufficient analyticity to derive the 
Regge theory — if anything, more so when Q 2 < 0. The analyticity properties are determined by 
the masses of the physical hadrons, and it is crucial that none of these has zero mass. The fact that 
gluons perhaps have zero mass is irrelevant because they are not physical particles, even though often 
we pretend they are in calculations. 
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With no zero-mass particle exchange, there is a finite "Lehmann ellipse" in the t-channel, within which 
the t-channel partial wave series converges for t > Am 2 , where m is the proton mass, and for suitable 
values of s — or, equivalently, v if we are considering deep inelastic scattering. For example, for the 
amplitude T\(v, t, Q 2 ) whose imaginary part at t = is F\(x, Q 2 ) with x = Q 2 /2u, 

oo 

T l {v,t,Q 2 ) = Y J (^+l)ai(t,Q 2 )P l (co S 6 t ) (1.1a) 

t=o 

where v must be such that 

v-H 

cos0 t = 4 (1.16) 



^(lt + Q 2 )(\t-m2) 



has absolute magnitude no greater than 1. The Froissart-Gribov technique' 3 ] then uses a fixed-i 
dispersion relation to define an analytic function a(£,t,Q 2 ) which is equal to ai(t,Q 2 ) when the 
angular momentum £ takes the values 0,1,2,3,. . . , but is defined also for complex values of I and has 
suitable behaviour when \£\ — ► oo to allow a Watson-Sommerfeld transform of the partial-wave series 
(1.1a). This expresses the series as an integral*: 

with Pg(cos9t) the analytic continuation of the Legendre polynomial to complex i. The contour C is 
initially that of figure la, wrapped arround the positive real-£ axis on which are located the zeros of 
the denominator sm(ir£). But the properties of a(£,t,Q 2 ) allow it to be distorted to become parallel 
to the imaginary-^ axis, as in figure lb, so extending the range of values of v and t for which the 
representation of Ti(v,t,Q 2 ) is valid. 

From what is known about the analyticity structure, it is rather certain that the positions of the 
singularities of a(£,t,Q 2 ) in the complex £ plane do not depend on Q 2 , only on t. If we analytically 
continue in t to the value we need, t = 0, a singularity may try to cross the contour C, and we must 
distort the contour again so as to avoid this happening. Depending on whether the singularity is a 
branch point or a pole, we then have either figure 2a or 2b. In the case of a pole 



the integration around it yields a contribution 

ir[3(Q 2 ,t)P a(t) (cos9 t ) 



(1.4a) 



sin7ra(t) 

to Ti(u,t,Q 2 ). Now make an analytic continuation so that v becomes large, 

v > Qm (1.46) 

More correctly, the Froissart-Gribov technique splits the summation (1.1a) into two, one over even 
£ and the other over odd £. Each is converted into an integral, with an "even-signature" amplitude 
a + (£, t, Q 2 ) that equals a^(i, Q 2 ) for £ = 0, 2, 4, . . . and an "odd-signature amplitude a~ (£, t, Q 2 ) that 
equals ae(t, Q 2 ) for £ = 1, 3, 5, — 



2 




Figure 1: Contour C in the complex t plane for the integral (1.2) 
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(a) 



(b) 



Figure 2: Distortion of the contour C caused by (a) a pole, or (b) a branch point of a(£, t, Q 2 ) crossing 
it. 

Then | cos0 t | S> 1. For large z, 

Pi(z)~z £ (1.4c) 

so the Watson-Sommerfeld transform (1.2) becomes essentially a Mellin transform^ and the "Regge 
pole" a(t) contributes at t = 

&i(Q 2 K (0) (1.4d) 

to the large-i^ behaviour of T\{v, 0, Q 2 ). Here, b\{Q 2 ) is a constant multiple of f3(Q 2 ,0). In the case 
of T?,{y, 0,Q 2 ) its definition includes a kinematic factor which reduces the power of v by one unit, so 
since v = Q 2 /x this gives 

F 2 (s,Q 2 )~/(QV~ a(0) (l-5a) 

Regge theory gives no information about the function f(Q 2 ), other than that they are analytic func- 
tions with singularities whose locations are known! 6 ]. The power (1 — a(0)) is independent of Q 2 . 

In the case where the singularity that crosses the contour C is a branch point at i = a c (t), dragging 
a branch cut with it as shown in figure 2b, the simple power of x in (1.5) is replaced with 

r a c (0) 

F 2 (x,Q 2 )~ / d£f(£, Q 2 )x l ~ e (1.56) 

which results in F 2 (x,Q 2 ) ~ x 1-0 -^ times an unknown function of Q 2 and log(l/x). Again, the 
power (1 — a c (0)) is independent of Q 2 . One knows, from unitarity® , that if there are poles in the 
complex t plane, there must also be branch points. On the principle that it is usually the best strategy 
to try the simplest possible assumption first, and because there is good evidence that it is correct in 
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Figure 3: The coefficient functions /o and /i extracted from data at each Q 2 ; the error bars are from 
MINUIT 



purely hadronic interactions^ , in reference 4 we tested the hypothesis that the contribution to the 
structure function F 2 at small x from branch points is much weaker than from poles. We made a fit 
of the form 

2 

F 2 {x,Q 2 ) ~ ^2 fi(Q 2 )x~ e ' (1.6a) 

i=0 

We fixed the values of two of the powers from our knowledge' 7 ] of hadron-hadron total cross-sections: 

ei = 0.08 ("soft pomeron" exchange) , 
e 2 = — 0.45 (p, u>, f, a exchange) 

and extracted the value of eo from the small- x data — we used all data for which x < 0.07 and Q 2 
ranged from to 2000 GeV 2 . This gave e ~ 0.4 and provided an excellent fit to the data over the 
whole range of Q 2 . 

In our fit, we constrained each of the three coefficent functions fi(Q 2 ) to vanish as (Q 2 ) 1+ei as Q 2 — > 0. 
Gauge invariance requires the total F 2 to vanish linearly with Q 2 at fixed u, and these constraints 
ensure that each of the three separate terms in (1.6a) have this property. From the data, we found 
that the "hard-pomeron" term is present already at small Q 2 , though its coefficient function is very 
small, but it grows approximately logarithmically for large Q 2 . The soft-pomeron coefficient function, 
on the other hand, rises quickly with Q 2 for small Q 2 , peaks at a value of Q 2 between 5 and 10 GeV 2 , 
and then becomes small at large Q 2 . See figure 3. The data do not impose much constraint on how 
the meson-exchange coefficent function behaves for large Q 2 . 

Regge theory gives no quantitative information about the coefficient functions fi(Q 2 ), but at large 
Q 2 they should evolve according to perturbative QCD. That is, provided they correspond to leading- 
twist contributions, their behaviour at large Q 2 should be governed by the DGLAP equation' 1 ]. The 
constraint from Regge theory is that the positions of singularities in the complex I plane should should 
not change with Q 2 , and that new singularities do not suddenly appear as Q 2 increases. In our fit to 
the dataW , we have shown that this is fully supported by experiment. It is our purpose in this paper 
to discuss it further. 
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2 The DGLAP equation: general theory 

We write the singlet DGLAP equation as^ 







dlogQ' 

where 

Introduce the Mellin transform 



-u(x,Q 2 ) = f dzTr(z,a s (Q 2 ))u(-,Q 2 ) (2.1a) 

J X z 



+ (2 . 16) 



so that 



U(N,Q 2 )= [ dxx N - 1 u(x,Q 2 ) (2.2a) 
Jo 

u(x,Q 2 )= j ^Lx- N XJ(N,Q 2 ) (2.2b) 

As we have explained, U(7V, Q 2 ) is essentially the complex-partial-wave amplitude a(£, t, Q 2 ) at t = 0, 
with I = N + 1, so that a power contribution 

/(Q 2 )x- e (2.3a) 

to F 2 (x, Q 2 ) corresponds to a pole 

in U(iV, Q 2 ). If the small-x behaviour of F 2 is (1.6a), the rightmost singularity of TJ(N, Q 2 ) in the 
complex A^-plane is a pole at N = eo, and the integration in (2.2b) initially should be along a line 

Re iV = iV with N > e (2.4) 

It is evident from its definition (2.2a) that U(N, Q 2 ) is analytic to the right of this line in the complex 
iV-plane, and for Re N —> +oo it goes to for physical values of x, that is < x < 1. In what follows 
we shall analytically continue U(iV, Q 2 ) to values of N such that Re N < eo; then it is no longer given 
by the representation (2.2a), which then diverges. 

Insert the Mellin-transform representation (2.2b) into the DGLAP equation (2.1a): 

j d ^-"lii^ V{N ' Q2)= l WQ>))U(1V,Q>) (2.5a) 

where 

n(N,x,a s (Q 2 ))= f dzz N ir(z,a s (Q 2 )) (2.56) 

If the rightmost singularity of U(iV, Q 2 ) in the complex iV-plane is a pole at N = e , with residue 
fo(Q 2 )) we may move the contour of N integration to the left past this pole and so pick up a contribution 
to each side of (2.5a) that behaves as x~ e ° . This gives the differential equation 

* V f (Q 2 )=n(6 ,as(Q 2 ))fo(Q 2 ) (2.6a) 



dlogQ 



with II(eo, cts{Q 2 )) the analytic continuation to N = eo of 

n(N,a s (Q 2 ))= [ dzz N *{z,a s {Q 2 )) (2.66) 
J o 

To leading order in as(Q 2 ), the difference between n(eo, cts{Q 2 )) and n(e , x, as(Q 2 )) behaves as 
x e ° at x = and therefore is negligibly small. 

A differential equation such as (2.6) actually applies to the coefficient function of each power term 
f(Q 2 )x~ e in u(x,Q 2 ), not just the leading one. That is, a term that behaves as any power at some 
value of Q 2 remains a power. However, there are apparent complications, arising from the fact that, 
to lowest order in as(Q 2 ), Tl(N,x,as(Q 2 )) is singular at N = 0. We address this matter in the next 
section. 

In lowest order, the splitting function n(z, as(Q 2 )) introduced in (2.1a) is simply proportional to 
ots(Q 2 ), and likewise its Mellin transform tl(N, as(Q 2 )) of (2.6b), so that (2.6a) reads 

r^f (Q 2 ) = B(e )fo(Q 2 ) r = log |J (2.7) 

The solution to (2.7) is 

f (Q 2 ) = exp (T(Q 2 )B( eo )) fc(Q§) T(Q 2 ) = log (2.8) 

In terms of the standard anomalous-dimension matrix 7(j), 

B(N) _ c (% q (N-l) 2N f%G (N-l)\ 6 

The elements of j(j — 1) are plotted in figure 4.9 of reference 8 whose authors have kindly provided 
us with the output of their computations * , from which in figure 4 we plot the elements of the matrix 
B(7V) in the case Nf = 4, and in figure 5 its eigenvalues X±(N). 

If we expand fo(Q 2 ) as a linear combination of the two eigenvectors of B(e ): 

f (Q 2 ) = o +(Q 2 )v + (e o ) + <P -(Q>-(eo) (2.10a) 

where 

(K ± (N))- 2 = (B, o (A0) 2 + (A±(JV) - B m (N)) 2 (2Mb) 

then 

, 2 \ A ±( e o) 



50± 



2 f \ogQ 2 /A 2 Y^' 2 

{q) = Wq 2 /a 2 ) MQo) (2 - 10c) 

According to figure 5, one eigenvalue A_ is negative, so for large Q 2 only the first term on the right- 
hand-side of (2.10a) survives, and then the hard-pomeron coefficient functions satisfy 

f 0q (Q 2 ) = B qG (e ) 

f 0G (Q 2 ) A+(e ) - B qq (e ) { " ' 



* It is not made clear in the book that the qg plot includes the necessary factor 2Nf. 
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Figure 5: Eigenvalues X±(N) of the matrix B(iV) 
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when Q 2 is large enough. 

In our previous fit to the structure-function data we used the forms 



^ 2 )^°(e^) 1+, °( 1+xlog ( 1 + S)) 

/ o 2 \ 1+ei 1 

where the soft-pomeron power t\ was constrained to be 0.08 and the meson-exchange power e 2 was 
fixed at -0.45. This resulted in a \ 2 just l ess than 1.02 per data point for the 539 data points for which 
x < 0.07, including the real-photon total cross-section measurements for which y/s > 6 GeV. The fit 
gave the "hard-pomeron" power eo to be about 0.4 for which, according to figure 5a, A+(e ) is close to 
2. On the other hand, figure 5a would give A+(ei) somewhat larger, which is in obvious conflict with 
the behaviour of fi(Q 2 ) at large Q 2 extracted from the data: see figure 3b. This reinforces our view 
that the soft-pomeron term is higher twist. In a new fit, we have therefore left its form unchanged, 
and replaced the hard-pomeron coefficent function with 

«« 2 » = ^(q^) 1+< ° ( 1 + M 1 + ^f T < 2 - 12t > 

so there is one fewer parameter. In principle, if the meson-exchange term is leading twist (which 
presumably it is) we should also include a logarithmic factor in its high-Q 2 behaviour, raised to power 
A+(e 2 ) obtained by analytically continuing A+(iV) to negative N. However, this term gives only a 
very small contribution to the high-Q 2 data points and so we have left it as in (2.12a). We find that 
making the change (2.12b) increases the \ 2 ver y slightly, to a little over 1.05 per data point. 

Our fit omits any contribution from a singularity at N = 0. We now explain why this is justified. 



3 Singularity at N = 

To lowest-order in as the splitting function ir(z, as(Q 2 )) nas a z~ l singularity at z = 0, so that its 
Mellin transform Tl(N, cts(Q 2 )), defined in (2.5b), has an iV -1 singularity at N = 0. This means that 
every power x~ e , whether e is whether positive or negative, generates through (2.5a) an additional 
singularity in XJ(N, Q 2 ) at N = as Q 2 increases. The evolution of this N = singularity is not 
controlled just by the small-x behaviour of the structure function. 

To investigate these matters, write the lowest-order splitting function as 

w(z, a s (Q 2 )) = [p(z) + Po^" 1 ] (3.1) 

where p(z) is regular at z = 0. Inserting this into (2.5) gives 

/ ^-"m^^^^j ^-{P(iv,,) + ^P„(1-X«)} UW Q 2 ) (3,) 
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where P(N,x) is a Mellin transform of p(z) like (2.5b). It vanishes so rapidly at x = that it is 
a good approximation to replace it with P(A, 0). Because U(A, Q 2 ) is analytic to the right of the 
contour of the iV integration and, from its definition (2.2a), vanishes when |A| — ► oo, the last term 
inside the curly bracket integrates to zero, and the inverse Mellin transform of (3.2) is just 

U(A,Q 2 ) ~a 5 (Q 2 ) P(N,0)) + N^po U(A,Q 2 ) 



dlogQ 



us(Q 
2vr 



2\ 



P(A)U(A,Q 2 ) (3.3) 



Evidently a pole (2.3b) at A = e in U(A,Q 2 ) gives a pole in (d /d log Q 2 )V(N,Q 2 ). By taking the 
residue at the pole of both sides of (3.3) we see that the coefficient function f (Q 2 ) satisfies a differential 
equation like (2.6a), involving the analytic continuation in A of P(A) to A = e. This is true for either 
sign of e. 



The solution to (3.3) is 



U(A,Q^)~exp Clogg||^P(iV)]u(iV,Q 2 ) (3.4) 



where the constant C is defined in (2.9). It is evident that, to use this equation, we must make some 
assumption about XJ(N,Qq) at some value Qq of Q 2 . The common assumption^ is that, for suitable 
choice of Qq, it is regular at A = 0. Indeed, there is no reason to suppose that, for Q 2 chosen to be 
or very small, XJ(N,Qq) is singular at A = 0. That is, although the real-photon Regge amplitude 
a(£, t, Q 2 ) has a wrong-signature fixed pole' 10 ' at I = 1, there is no reason to believe that at small Q 2 
it has also a right-signature singularity in the spin-averaged amplitude. Nevertheless, as soon as Q 2 
becomes large enough for the evolution to be correctly described by the DGLAP equation, we learn 
from (3.4) that U(N,Qq) acquires an essential singularity at A = 0. A function such as U(A, Qq) 
that is analytic at = for a finite range of Q 2 cannot suddenly acquire a fixed singularity at N = 
when it is analytically continued in Q 2 . 

The explanation* of this contradiction is that it is wrong to use the lowest-order approximation to 
the splitting function, II (iV, x, as(Q 2 )), with its A -1 singularity, for any purpose where its behaviour 
near A = is important. The higher-order terms in its expansion involve more and more powers of 
ag/A and near A = these powers become larger and larger, so that the expansion cannot possibly 
converge. Our general argument, that new singularities cannot suddenly appear in U(A, Q 2 ) as Q 2 
increases, strongly suggests that, if one were to resum the powers of as/N one would find that in fact 
U(N,x,a s (Q 2 )) is regular at A = 0. 

This is suppported by what is known about this resummationl 11 ] . The DGLAP anomalous dimension 
7gg(A — 1) in (2.9) is equal to the BFKL anomalous dimension 7(A) calculated to all orders in as/N. 
It is obtained by solving the equation 



where x( x ) is the Lipatov characteristic function. To lowest order in as 

x(x) = —Zje — ip{x) — — x) (3.56) 



* We are grateful to Guido Altarelli and Stefano Catani for helping us to understand this 
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Figure 6: Plot of the lowest-order approximation (3.5b) to the Lipatov characteristic function 

which is plotted in figure 6. From this plot, one sees that there are two complex-conjugate solutions 
for j(N) at N = 0; these are finite. If one expands the solution to the implicit equations (3.5) for 
j(N) in powers of as, the terms in the expansion are singular at N = because the digamma function 
Y>(7) is singular at 7 = 0, but this expansion is completely invalid near N = 0. 

We now know! 12 ! that the lowest-order approximation (3.5b) to the Lipatov characteristic function is 
inadequate. When the next term in the expansion is included the curve in figure 6 is pulled downwards 
and j(N) is no longer complex at N = 0, but it is still nonsingular, as our general argument has 
demanded. However, because the second term in the expansion of the Lipatov characteristic function 
is so large, it is likely that subsequent terms are too, and so we do not yet have a reliable estimate of 
just how large j(N) is at N = 0. But we know that it is nonsingular. 

Therefore, that in fitting data one should not include in the splitting function any singularity at N = 0. 
The only iV-plane singularities are those already present in the structure function at small Q 2 . These 
are the standard singularities of Regge theory — the soft pomeron, the mesons, and possibly also a 
hard pomeron. There may well be branch-point singularities, but the data are consistent with the 
assumption that their contribution is relatively small at presently accessible values of x and Q 2 . This 
is the approach we followed in our previous paper' 4 ', and we found that it was highly successful. We 
found also that the usual attitude to higher twists is more than a little suspect: rather than making 
a fit that starts by assuming they are as small as the data will allow, one should recognise that the 
data strongly suggest that a significant part of the structure function at Q 2 values of 5 GeV 2 or less 
is higher twist — the soft-pomeron term, and maybe also meson exchange. 

This research is supported in part by the EU Programme "Training and Mobility of Researchers", Networks "Hadronic 
Physics with High Energy Electromagnetic Probes" (contract FMRX-CT96-0008) and "Quantum Chromodynamics and 
the Deep Structure of Elementary Particles" (contract FMRX-CT98-0194) , and by PPARC 
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